Theory of which path dephasing in single electron interference due to trace in 

conductive environment 



Pawel Machnikowski 1 '^ 

institute of Physics, Wroclaw University of Technology, 50-370 Wroclaw, Poland 

A single-electron two-path interference (Young) experiment is considered theoretically. The de- 
coherence of an electron wave packet due to the which path trace left in the conducting (metallic) 
plate placed under the electron trajectories is calculated using the many-body quantum description 
of the electron gas reservoir. 
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I. INTRODUCTION 

Understanding the interaction between a microscopic 
quantum system and its macroscopic environment is one 
of the most exciting challenges of the modern quantum 
mechanics. Such an interaction is known to perturb the 
phase relations between the components of quantum su- 
perpositions which damages the system ability to show 
quantum interference. This effect, known as decoherence 
or dephasing, is essential for our understanding of the 
origins of classical behavior. 1 '- Moreover, it limits the 
feasibility of exploiting the quantum properties of micro- 
and mesoscopic systems for performing useful tasks, out 
of which quantum computation 3,4 is perhaps the widest 
known. 

Dephasing phenomena manifest themselves in any 
quantum system. It is, however, natural to seek for 
their realization in the most clear form. Such an ul- 
timate demonstration seems to be a realization of a 
"thought" experiment opening many quantum mechanics 
textbooks^: vanishing of the single-particle interference 
fringes in a two-slit (Young) experiment due either to 
some perturbation of the environment state or to a con- 
trolled measurement. Such a trace (often referred to as 
which path information) left in the measurement device 
or in the environment allows one (at least in principle) 
to identify the path chosen by the particle and destructs 
the interference to the extent depending on the distin- 
guishability of the paths. 

Experiments of this kind were realized with electrons 
in a semiconductor systemSiL£ as well as in optical* 9 ^ 
neutron^i and atomioiSiiiii 4 * 1 ^ interference setups. It 
was also suggested— that an analogous demonstration 
in the time domain could be based on interference 
experimentsiiiS*^ in semiconductor quantum dots. 

Another realization has been proposed^fiiSi involving 
free electron interference with a plate of finite conduc- 
tivity placed below the electron paths and parallel to 
them. A model for the decoherence effect in such a setup 
has been formulated on phenomenological basis^ using 
the picture of image charges screening the external flying 
electron. The image charge would move parallel to the 
real electron inducing energy dissipation (Joule heat gen- 
eration) due to finite conductivity of the material. The 
dissipation effect was quantitatively estimated using an 
earlier calculation of velocity fields penetrating a metal— 
The essential qualitative features of the resulting dephas- 



ing effect is the dependence on the distance to the plate 
(which governs the overall intensity of the which path 
trace) and on the distance between the paths (on which 
the degree of distinguishability depends). 

With the state-of-the-art experimental techniques, 
such a free-electron interference and controlled de- 
phasing experiment is feasible. It is possible to ob- 
serve interference fringes in experiments with single 
electrons ^SliSk It has been proposed to use a single- 
electron interferometer— to show the controlled dephas- 
ing effect due to the conducting plate j 27 ' 28 

The present paper aims at the development of a fully 
quantum description of the electron dephasing under the 
specific experimental conditions as described above. On 
the microscopic level, the which path information is trans- 
ferred from the electron (the system) to the conducting 
plate (the environment, or reservoir) by exciting the elec- 
trons in the plate. This effect may be described by ana- 
lyzing the electron coupling to the electromagnetic modes 
affected by the presence of a conducting surface, as pre- 
viously done for a superconducting 29 and conducting 30 
plate. In order to capture the essential (dissipative) ef- 
fect of the excitations in the reservoir it is, however, prof- 
itable to formulate the description in terms of the direct 
coupling between these two systems which is achieved by 
using the Coulomb gauge. 

A description developed in a different context 31 ' 32 
shows that such a which path decoherence effect can be 
quantitatively represented as the overlap between the 
spectral density of the reservoir fluctuations and the ap- 
propriate spectral function related to the unperturbed 
evolution of the system. The present calculation follows 
the same path, with the reservoir fluctuations expressed 
in terms of the standard longitudinal dielectric function 
of the conductor. It turns out that the interaction be- 
tween the flying electron and the charges in the plate 
involves low frequencies but high momentum transfer, 
which suggests that description going beyond the Debye 
model is necessary. Therefore, the present calculation 
includes a complete model of the quantum properties of 
the reservoir. 

The paper is organized as follows. In the following 
Section [H] the system under discussion and its model are 
presented. Next, in Sec. lIIII the general framework of the 
theory is described. In Sec. lIVI the problem is reduced to 
finding the low frequency longitudinal dielectric function. 
This is done for a metallic plate in the following Sec. [V] 
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FIG. 1: Schematic plot of the system. The information about 
the path chosen by the electron is transferred to the plate by 
exciting the electron gas (local heating). 

where the quantitative results are found. The final Sec. 
IVII concludes the paper. 



II. THE MODEL 

In the proposed experimental setupf2L2£ non- 
relativistic electrons with energies in the range of 150 
eV to 3 keV are emitted by an electron gun. This 
corresponds to velocities v ~ 7 • 10 6 to 3 • 10 7 m/s (i.e., 
roughly one order of magnitude higher than the Fermi 
velocities for metals^). The electron beam is split into 
two paths separated by D ~ 10 -j- 300 /im. Both paths 
pass over a conducting plate of length L ~ 1 cm at 
a height z$ ~ 0.1 mm and then interfere on a screen. 
The part of the experimental system relevant for the 
dephasing effect is shown in Fig. ^ It will be assumed 
that over the plate the electron paths are parallel. The 
plate is taken to have a length L and a very large height 
H and width L y . The coordinate system is oriented in 
such a way that the z = plane coincides with the plate 
surface and the electron paths are parallel to the x axis. 
In order to refer the present calculations directly to 
the experiment, a single electron wave function will be 
chosen in the form of a Gaussian quantum wave packet 
travelling over the plate at a speed v in a superposition 
of the two paths. We will assume that the electron wave 
packets corresponding to the two alternative paths do 
not overlap and that their extension is much smaller 
than any distance in the experiment geometry. For the 
electron velocities as given above, the dispersion of a 
minimum-uncertainty electron wave packet with initial 
width of a few micrometers leads to an additional spread 
of ~ 0.01 pm per 1 cm of flight path, so the effects of 
dispersion may be neglected. It will also be assumed 
that the time interval between the consecutive electrons 
is much larger than the relaxation times relevant to the 
conducting plate so that each flight event is independent 
of the previous ones. The experiment is done at room 
temperature. 

As mentioned above, contrary to the previous 
work(22i2H ^ ne p re sent analysis uses the Coulomb gauge. 
In this way, the longitudinal electromagnetic field is elim- 
inated in favor of direct and instantaneous coupling to 



the electronic degrees of freedom in the metal. Thus, 
the problem is of longitudinal screening type and may 
be treated using the well-established knowledge on the 
density-fluctuation spectra of the electron gas and its 
relation to the longitudinal dielectric functioni^i This 
choice relates also the present calculation more closely 
to the original phenomenological concept (although the 
decoherence mechanism is described in a different way). 

The Hamiltonian for the Coulomb interaction of the 
electron in the beam with the electrons and ions in the 
solid is 

v = ^[ d ^[d^efMrl, (i) 

47re J J \r - r'\ 

where Pb{t) is the density operator for the electron beam, 
p(r) is the total charge density at a point r in the con- 
ducting plate (including both electrons and ions), e is the 
elementary charge, eo is the vacuum dielectric constant 
(the SI system of units is used). The ordering of opera- 
tors in Eq. Q results from the assumption that the elec- 
tron above the plate is distinguishable from those in the 
solid. Since a uniformly moving electron does not radiate 
transverse fields on macroscopic distances, we disregard 
the coupling to the electromagnetic vacuum via vector 
potential and calculate the contribution to decoherence 
only from the longitudinal excitations, thus neglecting 
the effect of the conducting surface via modification of 
the transverse electromagnetic field fluctuations ,SL22i It 
should be noticed that we aim at a full microscopic de- 
scription of the interaction, hence no image potentials are 
inserted "by hand" in the Hamiltonian QJ. The interac- 
tion between the electron and the charges in the plate 
(classically described as the image potential) will result 
from the microscopic theory, while the mutual screening 
of electrons in the plate must be included in the correct 
description of the longitudinal response of the interacting 
electron system. 

Since we consider a flight of exactly one electron over 
the plate, the corresponding density operator may be 
written in the form p&(r) = |r)(r|, where \r) is the posi- 
tion eigenstate of the electron. In the interaction picture, 
the Coulomb part of the interaction Hamiltonian reads 

Vc(t) = J d\uMr){r\U {t)^- J d\>^, (2) 

where p(r,t) — U^(t)p(r)Uo(t) and Uo(t) is the evolution 
operator for noninteracting subsystems. 

The electron may travel along one of two paths, with 
the corresponding quantum states 

17 (*)|0(1)> - J d 3 W> 0( i)(r,i)|r), (3) 

where ipo,i( r i t) are the electron wave functions in the po- 
sition representation. Since the essential decoherence ef- 
fect is related to path distinguishability which affects only 
the relative phase between these two quantum states, 
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processes leading outside the two-dimensional subspace 
spanned by these states may be neglected. Projecting 
the electron states onto this subspace one gets 



(0(l)|V o (t)|0(l)) 



47T60 

(1(0)|^)|0(1)) = 



tfV|</> (i)M)r / d A r 



r — r' 



Attcq 



d 3 r^ m {r,t)^ 0{1) (r,t) / dVp_l^- = 0, 



j3 ,p(r',t) 



\r — r' 



under the assumption that the wave functions corre- 
sponding to the two basis states do not overlap. Hence, 
we have 

Vcit) = |0>(0|-^/d 3 r / dV|MM)| 2 ^l (4) 



UXil 



47T60 



dPr / <£V|V>i(r,t)|- 



P(r',t) 



III. 



THE WHICH PATH DEPHASING OF THE 
ELECTRON BEAM 



In this Section the general framework for the descrip- 
tion of the perturbation to the electron state (in the 2- 
dimensional subspace defined above) is described. The 
approach is based on the second order expansion of the 
evolution equation for the density matrix in the pres- 
ence of the instantaneous Coulomb interaction with the 
charges in the conducting plate. As a result one gets an 
expression for the visibility of interference fringes in a 
single-electron interference experiment in terms of a set 
of correlation constants depending on the electron evo- 
lution and on the reservoir properties. Qualitative dis- 
cussion based on the general form of these correlation 
constants is also presented here, while the quantitative 
analysis is performed in the following Sections. 

According to the standard theory^ the final reduced 
density matrix of the electron subsystem may be written 
as 

g(t) = U (t)g(t)uUt), 

where p(t) is the reduced density matrix in the interac- 
tion picture. In the second order approximation (valid as 
long as the overall perturbation effect is weak) the latter 
reads 



g(t) 



Qo 



dTTr R [V c (T),Q Q ® qr] 



(5) 



dT , Tr R [Vc(T),[V c (T , ),Qo®e R ]]. 



Here Ttr denotes the trace over the reservoir degrees of 
freedom (i.e., electrons in the plate) and it has been as- 
sumed that at the initial time the system state is separa- 
ble into the product go® gR, where go is the initial state 



of the electron and g R is the thermal equilibrium state 
of the electron gas in the conducting plate (this results 
from the long time interval assumption). 

Since the average fluctuation of the electron gas density 
in the plate vanishes in equilibrium, the leading order 
contribution to the dephasing is the second term. It may 
be written in the form 

Ag = -i[h t , g ] - ^{A, g } + ®[g ], (6) 

where the first contribution is a unitary correction and 
the two other ones describe dephasing ({•,•} denotes an 
anti-commutator). In the present case it may be shown 
that ht ~ I, so that this term does not contribute. Phys- 
ically, it contains the correction to the electron motion 
due to the interaction with the image charge induced in 
the metal. Due to the symmetry between the paths, these 
corrections are identical for both states and do not induce 
any nontrivial evolution within the restricted subspace. 
Their overall effect may be estimated by considering the 
attraction force e 2 / (16tt£oZq) acting during the time of 
L/v, which results in a negligible vertical shift of at most 
100 nanometers. 

The operator A and the superoperator <!> are 



a=J2 R*Ml *[e]= E 



R^m\g\m, 



i=0,l 



ij=0,l 



where the correlation constants are 



Ri 7 



(7) 



f dr f dr' f d 3 n f d 3 r[ f d 3 r 2 f d 3 r' 2 

Jtn Jtn J J J J 



|^(n,T)| 2 |^-(r 2 ,T')| 2 



r 2 



(p«,r-r')pK)). 



Note that the integration in Eq. is symmetric in the 
time variables, while the antisymmetric part has been 
separated into h t in Eq. © (see Ref. [3^ f° r technical de- 
tails). Obviously, Rij = R*^, in fact, both these constant 
are real and equal to each other due to the symmetry be- 
tween the beams which excludes any relative phase shifts. 

Although due to the finite size of the conducting plate 
(p{r[,T)p(r' 2 )) does not have the full translational sym- 
metry, in the transverse (y) direction it may depend only 
on j/2 — Vx- Hence, if the wave functions ipo,i( r ) differ 
only by a shift along y then i?oo — Ru- The evolution 
conserves the diagonal elements of the density matrix, 
while the off-diagonal ones change according to 



(0|Ag|l) = AiJ<0|ft,|l) 



AR = Ro 



Rn 



Neglecting the plate edge effects (homogeneous ap- 
proximation: the decoherence effects simply accumulates 
while the electron is flying over the plate) and in the 
absence of reservoir memory (Markovian limit) the sys- 
tem state is separable at any time, so that the above 
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formula holds for any time step. Then one expects that 
AR should be proportional to the path segment Al trav- 
eled by the electron and the quantity A -1 = —AR/Al 
becomes the dephasing rate per unit path length. In 
this case, starting from the equal superposition state 
go = IMM \ipo) = (|0) + |l))/\/2, one gets after the 
flight over the plate of length L 



0=i(|O><O| + |l)<l|) 



l e -L/X 



11 + 11 



The detection probability for an electron at point r is 
proportional to 

J(r,t) = (r\g(t)\r) = (r\U (t)g(t)ul(t)\r) 
= \{\Mr,t)\ 2 + \Mr,t)\ 2 



+e 



- L ' x [tfS(r,t)Vi(r,t)+H.c.]} 



where we used the position representation defined in 
Eq. 10 At a point where the intensity j{r) of the two 
beams is equal, i.e., ^1,2(1",*) = \/j{r)e^ 1 > 2 ^ r ' t ' , one has 



I(r,t) = 2j(r) 



1 



-L/X 



cos A0(r) 



where A<p(r) — 4>i(r,i) — 4> 2 (r,t). In the two-slit exper- 
iment this phase difference does not depend on time but 
varies from point to point as a result of the difference of 
the corresponding paths lengths, which leads to the in- 
terference picture (it is assumed that j(r) varies slowly in 
space) . The visibility of the interference fringes, defined 
in the standard way using the maximum and minimum 
values of I(r), is then 



-L/X 



Before presenting the detailed calculations it may be 
worthwhile to discuss the effect qualitatively. First, if the 
correlations are strongly local, 

R(r 1 ,r 2 ,t) = 

1 ' '■ , (P(r[,t)p(r' 2 )) 



(47re /i) 2 
6{n - r 2 )R(r u t) 



J J |n-ri| 



r 2 - r' 2 



then 

R 



13 — 1 dr I dr' / d 3 r\ 
t J t J 

2|„/. /„ _\|2 ; 



since different wave functions do not overlap. In this case, 
AR does not vanish and a dephasing effect appears. 

On the other hand, for infinitely long-range correla- 
tions, one can write 

R(r l ,r i ,t) = f{r 1 )f(r 2 )R{t), 
where f(r) — f(x, z) describes the dependence of the 
interaction on the electron position relative to the plate. 
Then 



Rt 



t pt 
dr / dr' 

to •* to 



xR(t-t') 



dxdz<fr(x, z 1 t)/(x, z) 
dxdz<&(x, z, t )/(#, z) 



where 

$(x,z,t) = J dy\ip (r,T)\ 2 = J dy\i/>i(r, r)| 2 , 

since the wave functions differ only by a shift along y. 
Therefore i?oi = R11 an d the dephasing effect vanishes. 

In order to find quantitatively the decoherence path 
A the correlation constants Rij must be evaluated using 
the correct description of the properties of the electron 
gas in the conducting plate. This is the subject of the 
following Sections. 



IV. THE CORRELATION CONSTANTS FOR 
THE ELECTRON GAS SCREENING RESPONSE 

The goal of the present Section is to express the cor- 
relation constants defined in Eq. Q in terms of a stan- 
dard material response function (dielectric function) and 
a spectral function pertaining to the unperturbed evolu- 
tion (free flight) of a single electron over the conducting 
plate. Since the dielectric function is expressed in the 
momentum space and frequency domain, Eq. J7J must 
be first Fourier-transformed. 

As discussed in Sec^J the electron states are described 
by Gaussian wave packets moving with the velocity v, lo- 
calized along two parallel paths, separated by a distance 
D, at a fixed distance Zq from the plate, 



1 



1 (x-vty , (v- Vi r (z-zqV 
"2 — ii — + — m — + — m — 



for i = 0, 1, where j/0,1 = ±-D/2 and l x ,y,z are the wave 
function widths in the three directions. For such a Gaus- 
sian state one has 



/ 



d 3 r\Ur,t)\ 2 



I 

2^2 



d 3 p 



-ipr e i(pxVt+p y yi+PzZ ) g- 
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The conducting plate occupies the volume defined by — L/2 < x < L/2, —L y /2 < y < L y /2, —H < z < 0, where L 
is the finite plate length while L y and H are assumed to be very large. Using the fluctuation-dissipation theore m 36 ! 37 , 
the density-density correlation function may be expressed by the imaginary part of the full dielectric function of the 
conducting plate 



flip 1 \ "* 2 

ire 2 V ^ Q 



-^coth^Im 



1 



e{q,oS) 



(9) 



where the expansion in space is made in terms of the discrete set of plane waves in the finite volume of the plate (with 
periodic boundary conditions). 

Let us assume that the total dephasing accumulated during the whole flight over the plate is weak. In this case, 
the perturbative Eq. JSJ yields a good approximate description of the effect. Since the interaction takes place only 
when the electron is over the plate, one may shift the initial and final times of the evolution to and t in Eqs. J5j and 
(0 to — oo and +00, respectively. Then, substituting Eqs. JHJl and into Eq. JjJ and performing the integral over 
times we get 



R; 



167r 5 fte V 



duj coth — Im 



e(q,(J) 



d z r' 



d 3 p f d 3 p' 



x6(p' x v — w )e i [ _ ( p ' _<3 ^ , '2+PHW+P^ Zo ]e _ 3[(^Pi) 2 +( Z BPB) 2 +(^Pi) 2 ] . 

Since in the y direction the plate is very long, the integrals over y[ 2 yield a Dirac S. On the other hand, in the x 
direction, the plate length is limited. Performing the four integrations over x' x 2 , y[ 21 followed by those over p X:V and 
p' x one arrives at 



RtJ 4ir 3 e hv 2 V 



duj coth — Im 



1 



£{q,Lu) 



(10) 



xe 



'.qy{Vi-Vj) e -l[{^T+ih<lvfy 



2wLS L (q x + -)\L 



Here <5l(<?) 



and 



L = 



dp z 



1 



{uj/v) 2 + ql+p\ 



-e 4 



dze~ i ^ Pz+q ^ z+iPzZa 



-H 



H—>c 



q2q 



where q = [q 2 + (uj/v) 2 ] 1 ^ 2 . In this result, terms pro- 
portional to (l z /z ) exp[—(zo/l z ) 2 ] have been neglected, 
since l z -C zq. 

Since the characteristic momentum scales of the system 
are at least of the order of 1/zq and z n -C i, the broad- 
ening of the function Sz(q), which is of the order of l/L, 
may be neglected and one can write 8l{q) ~ 8(q). In this 
way one neglects the corrections related to the approach 
to the edge of the finite plate and to the fly-away phase 
after crossing it, compared to the dephasing accumulated 
during the flight directly over the plate. It should be 
noted that Eq. IjlUfl is limited to weak dephasing but 
it does not involve any Markovian approximations. On 
the other hand, the memory time of the screening re- 
sponse of the electrons in metals is of order of inverse 
Fermi energy, i.e., femtoseconds, which is many orders of 



magnitude shorter than any time scale of the problem. 
Vanishing memory together with the proportionality of 
dephasing to the length travelled over the plate allow us 
to interpret the perturbative result as a dephasing rate 
(per unit path length) and to obtain an exponential de- 
cay as a solution of the corresponding rate equation. As 
discussed in the Appendix^ the same result is obtained 
by coarse-graining the flight over the plate and explic- 
itly using the short memory assumption, which directly 
leads to correlation constants Rij proportional to time 
over each small time step. 



Using Eq. (|1(J|) in the above approximation and replac- 
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ing the summation over q by integration one may write 
e 2 L 1 



where 



AR 



duo hfiuo 
i — coth — — 
e hv (2tt) 3 J uo 2 



d 3 q 



Im 



1 



e(q,u) 



S(q,u), 



(11) 



S(q,co) 



5 sin 9 



qv 



[1 



^iDq cos <p cos #ol 



xe 



K^) 2 -^^?) 2 ™ 32 ^ 0082 S -[2z q+^(ql,) 2 y cos? <p cos? 8 + (^) 2 



cos 2 tp cos 2 6»o+ (^)' 



1 — Erf [ < / cos 2 ip cos 2 



, 2 



v 2g2 



where we write q = q (sin 9, cos 9 cos (p, cos 9 simp) and denote 9o = axcsin[uo/ (vq)]. Since zq is a large (macroscopic) 
distance, for q 3> uo/v ~ l/^o the above function is strongly peaked around </? = ±7r/2. On the other hand, the 
dielectric function of a conductor extends over momenta of the order of the Fermi momentum, <;~I;f> 1/zq and its 
imaginary part vanishes for low momenta, so that for the relevant momentum values one can write 



S(q, ui) = \s(q, u)8 (sin 9 + [s - |) + 6 (<p + |) 



(12) 



where 



S{q,u) 



27T rir/2 

/ cos9d9S(q,uo) 

ii J-tt/2 



Physically, the above formulas mean that the momentum transfer from the flying electron to the electron gas excitations 
in the (x, y) plane may be at most of order of 1/zq and the energy transfer of order of Hv/zq. 

The integral over 9 is trivial while that over ip may be performed by a saddle point approximation around ip = ±7r/2, 
writing cos(</? ± tt/2) w ±tp which, upon further substitution <p = [(qv/uo) 2 — 1] _1 / 2 m and extending the integration 
limits, leads to 



S(q,uo) = 







-1/2 

POO 




1 - 








J — CO 





-2iS^v / TT? 3 ' 



xe 



i(^) 2 -K^) v -M^) 2 ^+^l 

4 



2v 



The analysis of the spectral function S(q, uo) (see Fig. |5J) 
shows that it decays exponentially for uo > v/zq while for 
a fixed uo it varies with q on an interval ~ 1/zq to reach 
a plateau for q 3> w/w, where it attains a constant value, 
dependent only on uo, 



S(q,w\ 



q^>uj/v 



S(uo). 



The description simplifies considerably if one assumes 
that the size of the electron wave packet is much smaller 
than the distance to the conducting plate, i.e., l x ,y,z <C 
zq. Note that keeping the finite wave packet size leads 
only to negligible quantitative corrections. In particu- 
lar, it does not introduce any momentum cut-off. In the 



approximation l x ,y,z 0, the asymptotic value may be 
written in the form 



S(w) 



du 



1 - e 



— l^r- u 



1 + u 2 



(13) 



In the low frequency range selected by the function 
S(ui), one has always Im£ _1 (q,w) ~ ui (this is a gen- 
eral property^ see also below). Moreover, hv/(knZo) ~ 
1 K, so that at room temperature one may write 
coth(hfiuo / 2) w 2k-sT / (fouo) for all relevant frequencies 
lo < v/zq. Therefore, the frequency integral in Eq. i|ll|) 
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FIG. 2: The form of the function S(q,u>) for D/z = 0.1 and 
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FIG. 3: The "geometrical" function 7(1). 



will always have the form 



00 , S(uj) 

dbJ 



where 



7(3) = r 



16 



0(z 4 ), 



' 1 — 



■In 



4 1 + n 2 



(14) 



(15) 



[using the asymptotic form given in Eq. Ijl3[l]. This func- 
tion depends only on the geometry of the system. More 
specifically, the only relevant parameter is the ratio D/zg- 
The function 7(1) is plotted in Fig. EI 

As for the momentum integral in Eq. (|llf> . it obviously 
depends on the specific properties of excitations in the 
electron gas and may be very different depending, for 
instance, on the relation between the thermal energy k^T 
and the Fermi energy, as well as on the conduction band 
structure. In any case, however, this integral may be 
partly performed using Eq. (|12[l. 



d 3 q 



Im 



1 



e(q,uj) 



x—5 I sin# + 
2 V qv 





1 


j dq z Im 






e(q,uj)_ 



where it is assumed again that qv ^> u. Hence, from 
Eq. E3J, 



AR = 



(16) 



k B Te 2 L (D\ 1 



dq z Im 



e(q,u) 



Since, at low frequencies, £2(9,^) <C e\(q,u>), one may 
write 



Im 



e{q,oS) 



£2(q,u) 
e\{q,u))' 



=q y =0 



Below, the results for metallic plates are presented. 



V. DECOHERENCE OVER A METALLIC 
PLATE 

In this Section, the decoherence effect for a metallic 
plate is considered. For the description of the reservoir 
response properties, we will use the quasi-free electron 
(isotropic) longitudinal dielectric function in the random 
phase approximation (the Lindhard theory) By assum- 
ing periodic boundary conditions this dielectric function 
may be applied to a finite system with negligible cor- 
rections as long as the typical momentum transfers are 
much larger than the inverse system size, which is the 
case in our problem. Obviously, the screening response 
originates only from the finite volume of the plate, which 
has already been included in the theory developed above. 

On the other hand, the Lindhard theory assumes that 
single-particle states are just plain waves confined in a 
box with infinite potential walls, which is not quite the 
case for a real piece of metal. In particular, the dielectric 
discontinuity at the surface acts as a trapping potential 
for carriers^ Since the surface states are confined in the 
vertical direction, low energy excitations must involve in- 
plane momentum transfer, which is not favorable in our 
problem. Thus, one should not expect large dissipation 
(dephasing) from these modes an their contribution will 
be dominated by virtual high-frequency excitations (sur- 
face plasmons) which will tend to screen the potential of 
the external charge, thus decreasing the dephasing. This 
effect will be neglected in the present analysis. 

Even in this quasi-homogeneous approximation, the 
strictly correct description of the screening response re- 
quires solving the complicated many-body problem, in- 
volving the electrons in the metal, the positive ions, and 
the mutual screening interaction between these subsys- 
tems. Under the conditions of the present study the 
treatment may be, however, greatly simplified with some 
loss of precision but retaining the essential features and 
yielding reasonable quantitative estimations. 

Thus, the calculation of the total dielectric function of 
a metal will be done within the approximation 37 consist- 
ing in simply adding the susceptibilities of the electron 
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and ion subsystems, 

e(q,uj) = l + XeM, w) + x P h(q, w) 



[i + x P h(g,^)] 



l 



1 + X P h(<7,^) 



Xei(q,w) 



where x c i(q, w), Xph(<Z, ^) are the electron and phonon 
electric susceptibilities, respectively. 

The phonon susceptibility, including the screening of 
the ion interaction by electrons, is given by2£ 



pi 



u*(q) 



where uj v \ = Ze[n/ '(eo-M)] 1 / 2 is the ion plasma frequency 
(here Z is the ion charge, n is the ion concentration and 
M is the ion mass). Since there is no momentum cut-off 
in the problem under discussion, the system response is 
dominated by excitations with q of order of the Debye 
wave vector kp), so that we may replace to(q) by its short 
wavelength value for which the standard simple theory 
yields Lu p i. Also, since u> ~ v/zq <C Dpi, we may put u> — > 
in the denominator. In this approximation, the effect 
of screening by the lattice excitations is fully contained 
in the ion dielectric constant 

ei = l + X P h(fc D ,0)w2. 

In this discussion the imaginary part of the lattice re- 
sponse describing the dissipative effects has been ne- 
glected. Such effects must always be resonant, i.e., they 
involve phonons with low frequencies to ~ v/zq. The 
density of states for such long-wavelength phonons is 
very low. Moreover, at the corresponding low momenta 
v/(c s zo) <C kp, where c s is the sound speed, the real part 
of the total dielectric function is very large due to electron 
contribution (see below), so that such low- frequency lat- 
tice excitations are additionally very strongly screened. 
Therefore, the lattice contribution to the dissipative pro- 
cesses [i.e., to the £2(5, w) function] is negligible. In terms 
of the general picture of the dephasing due to the trace 
in the environment this means that most of the phonons 
have frequencies much higher than that characteristic of 
the external field so that they follow the perturbation 
adiabatically, returning to the original state after the fly- 
ing electron is away and thus registering no trace of its 
passage. 

The total dielectric function including electron transi- 
tions within the conduction band and the lattice screen- 
ing as described above is therefore described by the for- 
mula 



s(q,u) = £i 



1 

1 + -Xci(<?,w) 



(17) 



where we take into account that in the isotropic model 
the dielectric function may depend only on the value but 
not on the direction of q. The electron susceptibility in 
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FIG. 4: (a) Solid line: the "material" function fi(x) for met- 
als with the values corresponding to a few metals (assuming 
ej = 2); dashed and dotted lines: first and second order ap- 
proximation according to Eq. 12411 . respectively; points: the 
values for the five metals corrected for the Coulomb and ex- 
change correlations within the Hubbard model, (b) The in- 
verse decoherence length as a function of the D/zo ratio for 
gold at T = 293 K, for three values of the electron energy 



the quasi-free electron model valid for arbitrary q is given 
by the Lindhard formula^ 



47r 3 eoq 2 



d 3 p 



n F (E p ) - n F (E p+q ) 
E p - E p+q + (j + i0+ 



(18) 



where E p — p 2 /(2m), m is the electron mass and n F (E) is 
the Fermi-Dirac distribution. Since the function S(q,uj) 
has an exponential spectral cut-off at frequencies much 
lower than the Fermi energy, one needs only the low- 
frequency part of the dielectric function (as already re- 
marked in Sec. IIV|) . To be specific, hv/z ~ lO -5 ^, 
where E-p is the Fermi energy (the relevant frequency 
range corresponds to radio frequencies). Thus, for the 
imaginary part of the susceptibility [Eq. l|TSJ|] we write, 
using the standard results^ 



t c<j for uj/v-p < q < 
otherwise. 



(19) 



The electron velocities used in the experiment satisfy v 3> 
Up for practically all metals. Thus, the condition q 3> 
lu/v is automatically satisfied and the asymptotic form 
of the spectral function given by Eq. (|13|) may always be 
used. 

The leading low-frequency term of the real part of 
Eq. CHI) is 



Rexei(<7,0) 



me 2 kp 

2TT 2 € fl 2 q 2 



1 



4fc| - q 2 
Aqk F 



In 



q- 



2kp 



2kv 



(20) 

where kp is the Fermi wave vector. 

The momentum integral in Eq. (|16f) may be written, 
using the explicit expression l|20|) . in the form 



/ 



dq. 



£?(«*, 0) 



w ~3, V 1 



2'KtQt\W , kp 



The lower value of the integral over q may be extended 
from u>/v to 0. This is a negligible correction since lo/v <C 
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FIG. 5: (a) The visibility of interference fringes as a func- 
tion of the path distance from a gold plate of length L = 1 
cm, assuming the separation between the paths D = 10 /im, 
for three values of the electron energy, (b) A simulation of 
interference fringes in an experiment with an electron beam 
spread along z, in which the fringes for a range of zo values 
are obtained simultaneously. Here D = 10 fira, E = 150 eV, 
xs is the coordinate along the screen. 



kp and the expression under the integral is regular at 
q — » 0. The function 



n{x) 



x 

T 



du 



Ami 2 



u\ 1 
In — 



2u 



(21) 

-i -2 



describes the properties of the metallic reservoir (the ma- 
terial properties of the system). It is depicted in Fig.QJi. 
Inserting the above result into Eq. I|16|) and using Eq. i|14|) 
one gets the characteristic dephasing path length for a 
metallic reservoir 



A~ 



AR 
L 



7 



(22) 



This result is shown for gold in Fig. for a few elec- 
tron energies. The dependence on material parameters 
is contained in the function fi which, for many metals, 
differs only slightly from that corresponding to gold (see 
Fig-Et)- For electrons with energy 150 eV the decoher- 
ence effect should be noticeable with a 1 cm plate already 
for D/zo ~ 0.1. In Fig. |5j the visibility of interference 
fringes for a specific system setup is shown, along with a 
simulation of the fringes. 

A standard way to quantitatively improve the descrip- 
tion of the screening response of an electron gas is to in- 
clude correlations between electrons leading to Coulomb 
and exchange hole around an electron. This can be done 
within the Hubbard model^ where the electron suscep- 
tibility (|18fl is replaced by 



(H) / > Xel(g^) 



G(g)xei(g,w)' 



(23) 



with 37 



1 



G(q) = 7T 



2q 



where A;tf is the Thomas-Fermi wave vector for a given 
metal. It turns out that this correction affects the value 



of fj, for typical metals only to a very little degree, as 
shown in Fig. 01a) (points). 

It may be interesting to obtain a closed formula for the 
decoherence length A or for the decoherence time Td and 
compare it to the phenomenological prediction of Ref.l20l 
As can be seen in Fig. for most metals on has 



me 



1. 



27reo£ifi 2 fcF 

In this range of values the asymptotic expansion is valid 

1 



fj,(x) w -x- —x- 



(24) 



(see Fig.0^). Retaining only the linear term one finds 



dq. 



e 2 {q z ,u) _ 7r m 



z =-2 



e\{q z ,Q) Ac-Mf 



which leads to the approximate formula for the decoher- 
ence length A, 



, k B T me 2 (D 
A = 'y — 



(25) 



For D < zq we use the expansion from Eq. (|15J) and write 
for the decoherence time Td = X/v 



32 r VA 



D 



(26) 



where AdB = 2irh/ \J2m,k B T is the thermal de Broglie 
wave length of an electron and 



(27) 



27reoe;ZQ/ifcF 



The formula 
given in Ref. 



2611 is essentially of the same form as that 
2CI However, here AdB refers to the elec- 
trons in the metal and involves the mass of these elec- 
trons. In fact, decoherence results from the perturbation 
of the reservoir state and not that of the systemi. Since 
the reservoir responds to the electromagnetic field of the 
flying electron, its state may depend on its charge but 
not on the mass. Therefore, it should be expected that 
particles of different mass but equal charge will undergo 
the same decoherence. 

The time constant r r was originally interpreted as the 
energy dissipation rate, i.e. the rate at which energy is 
transferred from the electron to the reservoir. A simi- 
lar interpretation is possible here, within a semi-classical 
approach. The electrons in the metal involved in inter- 
actions with external fields are those near the Fermi sur- 
face. The relative rate at which their energy changes is 
therefore 

_ x _ 1 dE 
Tr ~ ~Ep ~dt' 
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The variation of the energy of electrons is due to their 
acceleration by the electric field £ of the flying electron 
the distance to which is roughly zo, so that 

dE . e£ e 2 

— — = mvv = mv — vp , 

at m 47reo£i^o 

where vp is the electron speed at the Fermi surface and 
the screening by ions has been included. Using the equal- 
ities Ep — h 2 k F /(2m) and vp = Hkp/m one arrives at the 
formula l(77l) . 

In spite of the formal similarity of Eq. (|26[) to the 
phenomenological decoherence rate obtained from the 
energy losses due to dc resistance and Joule heating 20 , 
the present theory actually describes a different mecha- 
nism of decoherence. The dc resistance-based description 
assumes implicitly that formation of the screening im- 
age charge is dissipationless (adiabatic), involving only 
virtual transitions, and therefore reversible. The irre- 
versible, dissipative processes take place only due to car- 
rier scattering as the image charge moves beneath the 
metal surface. Since, at high temperature, such scatter- 
ing is mostly due to phonons, the which path information 
is effectively stored in lattice excitations. In contrast, the 
quantum description presented here shows that already 
the process of formation of the image charge is to a large 
extent dissipative, even in the absence of carrier-phonon 
scattering. Since the reservoir response to the external 
charge involves only excitations of low frequency but of 
arbitrary momenta the appropriate dielectric function is 
different from the commonly used Drude limit. As a re- 
sult, there is no direct correlation between the decoher- 
ence effect and the dc conductivity or resistivity. Instead, 
for metals the theory predicts roughly inverse proportion- 
ality to the Fermi momentum which is the only material 
parameter entering the result (apart from the electron 
mass). Quantitative comparison shows that for noble 
metals the rate given by Eq. Q27[) is many orders of mag- 
nitude higher than that resulting from resistive dissipa- 
tion. Therefore, the Ohmic resistivity effect is of minor 
importance. 

VI. CONCLUSIONS 

The decoherence effect on a single electron travelling 
over a conducting plate was calculated in a model re- 
flecting the conditions of a currently performed experi- 
ment. The dissipative response of the electron gas in the 
plate, tending to screen the external charge, generates 
a trace in the plate. The resulting distinguishability of 
the electron paths destroys the electron's ability to in- 
terfere. In order to describe this effect quantitatively a 
fully quantum model of the mutual interaction between 
the electron and the conductive reservoir was formulated. 
The decoherence effect may be expressed by the spectral 
density of the reservoir fluctuations which, in turn, is 
related to its dissipative properties, i.e., to the imagi- 
nary part of the inverse longitudinal dielectric function. 



The resulting decoherence was described for a metallic 
reservoir but the qualitative form of the result depends 
only on the universal linear frequency dependence of its 
imaginary part. Therefore, the decoherence effect will 
be qualitatively similar for any plate material (e.g., for 
semiconductors) as long as the electron gas response is 
dominated by the low-frequency sector. Whether this is 
the case for a specific system, will depend on the interplay 
of the experimental conditions and material parameters. 

The results presented in the paper explain the mecha- 
nism of electron dephasing on a fully quantum level and 
give quantitative estimations of the effect, depending on 
the material parameters and the geometry of the exper- 
imental setup. Therefore, they should be helpful for ex- 
perimental studies in this fundamental field. 
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APPENDIX A: EXPLICIT MARKOVIAN 
APPROXIMATION 

In this Appendix, Eq. I|10[l is re-derived in a way that 
clearly displays the Markovian approximation by making 
an explicit use of the short memory time r mcm of the 
screening response which leads to correlation constants 
which are directly proportional to time. Now we consider 
the evolution during the time segment to < t < to + 
At during which the electron is moving over the plate. 
We assume that At is much longer than r mem but short 
enough for the accumulated perturbation to be small. 
Since T mem ~ h/ Ep is of order of femtoseconds and the 
time of flight over a 1 cm plate is in the nanosecond range 
such time-slicing is always possible in the experimentally 
interesting situation of non-complete dephasing. 

Let us start again by substituting Eqs. (JSJ) and @ 
into Eq. JJJ. Now, however, we first integrate over 
y'12, which yields 4ir 2 6(p y + q y )6(p' y — q y ), and over 

x{ 2 , which yields Air 2 S^\p x + q x )6^' ] {jp' x - q x ), where 

5^\q) = sm(qL/2)/(Trq). Based on the same argument 
as before, for L zq the latter may me replaced with 
a Dirac delta. The trivial integration over p x ,y and p' x y 
yields a formula containing the time integration in the 
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form and extended the limits of the second integration because 

t t Mqif) decays on time scales much shorter than At. 

dr dt ' e r i( i* v ( T - T '} 'M q (r — r') Using this result one obtains 

J tn J tn 



hi 



dr / due- tq * vu M q {u) 



2 

where we defined the memory function 



Rij 



2irAt coth — ^- Ime^g, "?*<o,? B ,? 



A/ ' J 1 £ ^coth^^Im 



u=\q x \v> 



poo j-o which is identical to the L — > oo limit of Eq. 1)10(1 with 

M q (t)= dwe-^coth-^-Ime'^q.uj) L = vAt. 

Jo 2 
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